High-resolution African HLA resource uncoversHLA-DRB1expression effects underlying vaccine response (preprint)

How human genetic variation contributes to vaccine immunogenicity and effectiveness is unclear, particularly in infants from Africa. We undertook genome-wide association analyses of eight vaccine antibody responses in 2,499 infants from three African countries and identified significant associations across the human leukocyte antigen (HLA) locus for five antigens spanning pertussis, diphtheria and hepatitis B vaccines. Using high-resolution HLA typing in 1,706 individuals from 11 African populations we constructed a continental imputation resource to fine-map signals of association across the class II HLA observing genetic variation explaining up to 10% of the observed variance in antibody responses. Using follicular helper T-cell assays, in silico binding, and immune cell eQTL datasets we find evidence of HLA-DRB1 expression correlating with serological response and inferred protection from pertussis following vaccination. This work improves our understanding of molecular mechanisms underlying HLA associations that should support vaccine design and development across Africa with wider global relevance. Teaser High-resolution typing of HLA diversity provides mechanistic insights into differential potency and inferred effectiveness of vaccines across Africa.

Estimating the relative proportions of SARS-CoV-2 strains from wastewater samples (preprint)

Wastewater surveillance has become essential for monitoring the spread of SARS-CoV-2. The quantification of SARS-CoV-2 RNA in wastewater correlates with the Covid-19 caseload in a community. However, estimating the proportions of different SARS-CoV-2 strains has remained technically difficult. We present a method for estimating the relative proportions of SARS-CoV-2 strains from wastewater samples. The method uses an initial step to remove unlikely strains, imputation of missing nucleotides using the global SARS-CoV-2 phylogeny, and an Expectation-Maximization (EM) algorithm for obtaining maximum likelihood estimates of the proportions of different strains in a sample. Using simulations with a reference database of >3 million SARS-CoV-2 genomes, we show that the estimated proportions accurately reflect the true proportions given sufficiently high sequencing depth and that the phylogenetic imputation is highly accurate and substantially improves the reference database.

Assessing uncertainty in the rooting of the SARS-CoV-2 phylogeny (preprint)

The rooting of the SARS-CoV-2 phylogeny is important for understanding the origin and early spread of the virus. Previously published phylogenies have used different rootings that do not always provide consistent results. We investigate several different strategies for rooting the SARS-CoV-2 tree and provide measures of statistical uncertainty for all methods. We show that methods based on the molecular clock tend to place the root in the B clade, while methods based on outgroup rooting tend to place the root in the A clade. The results from the two approaches are statistically incompatible, possibly as a consequence of deviations from a molecular clock or excess back-mutations. We also show that none of the methods provide strong statistical support for the placement of the root in any particular edge of the tree. Our results suggest that inferences on the origin and early spread of SARS-CoV-2 based on rooted trees should be interpreted with caution.

Synonymous mutations and the molecular evolution of SARS-Cov-2 origins (preprint)

Human severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is most closely related, by average genetic distance, to two coronaviruses isolated from bats, RaTG13 and RmYN02. However, there is a segment of high amino acid similarity between human SARS-CoV-2 and a pangolin isolated strain, GD410721, in the receptor binding domain (RBD) of the spike protein, a pattern that can be caused by either recombination or by convergent amino acid evolution driven by natural selection. We perform a detailed analysis of the synonymous divergence, which is less likely to be affected by selection than amino acid divergence, between human SARS-CoV-2 and related strains. We show that the synonymous divergence between the bat derived viruses and SARS-CoV-2 is larger than between GD410721 and SARS-CoV-2 in the RBD, providing strong additional support for the recombination hypothesis. However, the synonymous divergence between pangolin strain and SARS-CoV-2 is also relatively high, which is not consistent with a recent recombination between them, instead it suggests a recombination into RaTG13. We also find a 14-fold increase in the dN/dS ratio from the lineage leading to SARS-CoV-2 to the strains of the current pandemic, suggesting that the vast majority of non-synonymous mutations currently segregating within the human strains have a negative impact on viral fitness. Finally, we estimate that the time to the most recent common ancestor of SARS-CoV-2 and RaTG13 or RmYN02 based on synonymous divergence, is 51.71 years (95% C.I., 28.11-75.31) and 37.02 years (95% C.I., 18.19-55.85), respectively.